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Abstract 

It is shown that the Lagrangian reduction, in which solutions of equations of motion that do not 
involve time derivatives are used to eliminate variables, leads to results quite different from the 
standard Dirac treatment of the first order form of the Einstein-Hilbert action when the equations 
of motion correspond to the first class constraints. A form of the first order formulation of the 
Einstein-Hilbert action which is more suitable for the Dirac approach to constrained systems is 
presented. The Dirac and reduced approaches are compared and contrasted. This general dis- 
cussion is illustrated by a simple model in which all constraints and the gauge transformations 
which correspond to first class constraints are completely worked out using both methods in order 
to demonstrate explicitly their differences. These results show an inconsistency in the previous 
treatment of the first order Einstein-Hilbert action which is likely responsible for problems with 
its canonical quantization. 
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I. INTRODUCTION 



Canonical quantization is the oldest, most rigorous, non-perturbative approach to quan- 
tization. It demands no new hypotheses (which is especially important in quantum gravity 
because of lack of experimental guides) and rests completely on the classical general theory 
of relativity and conventional methods of quantum field theory. For a discussion of the 
problems that one faces in trying to establish a connection between classical gravity and 
models built on new hypotheses, see the review |1| . The first step in the canonical approach 
to any theory is to cast it into Hamiltonian form; analysis of this step is the main subject 
of this article. 

The search for a canonical formulation of the Einstein-Hilbert (EH) action began after 
initial developments in analyzing the dynamics of singular (gauge invariant) systems where 
constraints arise 0,0,0]. 

Almost immediately after Dirac presented his work^ on constraint dynamics the first 



attempt to apply his algorithm to t 



le gravitational field was made by Pirani, Schild and 



Skinner 0, , and by Dirac himself 

In the above mentioned articles, the metric formulation of the EH action was used 

5,(^"^) = J d'x^g^'-R,, (r, dT) , (1) 

where d is dimension of spacetime, g = det (gap), affine connections F^j, are equal to Christof- 
fel symbols = ^(?'*'°" {gij,a,u + gua,fi — gij.u,a) and R^i, is the Ricci tensor expressed in terms 
of F^j^ (see This is a "second order" formalism, as second derivatives of g^^, appear in 

Unlike 'ordinary' gauge theories, the Dirac analysis [9] cannot be applied directly to (0) 
because it is not known how to deal with second order derivatives using the Dirac procedure. 
(Both velocities and accelerations are explicitly present in (^.) To avoid this problem, the so- 
called gamma-gamma form L'^ was used as a starting point in obtaining the Hamiltonian 
for pure gravity 

L', (^?"^) = v^^7"^ (r^,F^^ - F^^F^,) . (2) 



^ The course of lectures_given at Canadian Mathematical Seminar, Vancouver, August-September 1949 and 
later published in 
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The Lagrangian of (j2I) differs from that in (^J) by a total divergence . The ehmination 
of such a term does not affect the field equations but the reduced Lagrangian of ^ is not 
relativistically invariant. (This was clearly stated in (5| and reflected in its title: "On the 
Quantization of Einstein's Field Equations" , not action.) The role of surface terms in the 
Hamiltonian formulation of General Relativity (GR) was discussed in |ll|. Recently, the 
peculiar features of surface terms were reconsidered from quite different perspective in 12 1 
where it was demonstrated that it is not possible to obtain the full EH action (only its 
gamma-gamma part) starting from the standard graviton action built from non-interacting, 
massless, spin-2 tensor field, and iterating in the coupling constant by having an interaction 
between the tensor field and its own energy-momentum tensor. 

The canonical approach based on instead of the EH action, is different from the 
canonical approach to ordinary gauge theories. If we use the invariance of original 
action is lost completely and not just its manifest form as in ordinary gauge theories after 
the time coordinate has been singled out. 

Dirac started his analysis using Q and later added the particular divergence term. Ac- 
cording to this is a logically incomplete procedure. He also introduced space-like sur- 
faces and fixation of coordinates in order to keep a space-like surface always space-like. This 
obviously destroys general covariance. In the conclusion to his paper which is, probably, 
not well-known, Dirac clearly stated what one gives up in his formulation: "One starts with 
ten degrees of freedom for each point in space, corresponding to the ten gfj_^, but one finds 
with the method here followed that some drop out, leaving only six, corresponding to six 
Qrs- This is a substantial simplification, but it can be obtained only at the expense of giving 
up four- dimensional symmetry. I am inclined to believe from this that four- dimensional 
symmetry is not a fundamental property of the physical world." In the next paragraph 
he continued: "The present paper shows that Hamiltonian methods, if expressed in their 
simplest form, force one to abandon the four- dimensional symmetry." (Italic of Dirac) This 
conclusion gives only the relationship between this simplest form of his Hamiltonian meth- 
ods and four-dimensional symmetry. Accepting Dirac's conclusion means that GR has to be 
finally reformulated without four-dimensional symmetry. This is what is done in Q] , where 
GR is reexpressed as a theory of evolving 3-dimensional conformal Riemannian geometries 
obtained by imposing two general principles: 1) time is derived from change; 2) motion and 
size are relative. 
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In contrast, if one believes that four-dimensional symmetry is a fundamental property 
of Nature and wants to keep this symmetry with the intention of eventually quantizing 
the EH action, one has to abandon the simplest Hamiltonian methods and try to find a 
Hamiltonian formulation that does not destroy four-dimensional symmetry right from outset. 
Two possible ways of doing this exist. The first one is to modify the Dirac procedure and 
work with the explicit dependence of the EH action on acceleration. The second is to find an 
equivalent formulation of the EH action that permits use of the standard Dirac procedure. 
(In addition to these orthodox approaches, there are a few more which are less developed; 
see p. 54 of and references therein.) 

In the first case, we can consider the EH Lagrangian as a Lagrangian with higher deriva- 
tives and try to apply the Ostrogradsky Hamiltonian formulation Q| with appropriate 
adjustments to accomodate singular systems. (It was clearly indicated by Ostrogradsky 
that he considered only non-singular cases.) The first systematic generalization to singular 
cases was given by Gitman and Tyutin (see also ^a|)- A full analysis of the EH action 
or some models where higher order derivatives enter only in such a way that they do not 
affect the equations of motion, to the best of our knowledge, does not exist. (The EH action 
in this respect is a kind of "one and a half" order system which probably creates problems 
in applying the Ostrogradsky method.) An attempt in this direction is due to Dutt and 
Dresden [l9|. 

The second approach which does not involve reduction of the EH action by the elimination 
of a total divergence makes the action first order in derivatives by introducing auxiliary fields. 
If by elimination of these fields, we can return exactly to the original action (including 
terms with second order derivatives), we have an equivalent form. This form of the EH 

n 

action is Einstein's affine- metric formulation J2CII|; it is just linear in first order derivatives, 
so the standard Dirac procedure can be applied similarly to the way it is applied to a first 
order formulation of ordinary gauge theories. Moreover, all terms of the first order action 
contribute to the equations of motion, as opposed, to the second order formulation, and so 
the effect of all terms can be studied on the same footing. 

Einstein considered g"^ and F^^ as independent fields without assuming F^^ = la/?}; 
since they are varied independently^. In j2^, he also proved that for symmetric (7"^ and F^^ 



This formulation was inspired by his search for unification of gravity and electromagnetism (he originaUy 
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this formulation is equivalent to (JJ and said that this was "the most simple and consistent 
way" of obtaining the field equations from the action principle. He also noted that with this 



formulation, no variation of fields on boundaries is needed (see also 22 L Appendix E). In 
this approach the division of variables into being bulk or surface, as in is avoided and 
all variables are treated on the same footing with field variations vanishing on the boundary. 
The Lagrange density of (eq.(3)) is given by 

where h"'^ = ^J—gg'^^ is just a simplifying notation and is not treated as an independent 
variable. (If we consider /i"^ = ^J—gg""^ as a change of variables, the functional Jacobian 
is field dependent in all dimensions d > 2 and for = 2 is singular.) Moreover, if we 



Sh° 



consider h'^^ as an independent field in (jH)) without taking into account the field dependence 
of the Jacobian, we cannot return to g'^^. 



lows from the solution of 
We then obtain the 



The equivalence of © to the second order form ^ {d > 2) fo 
the field equations for F^^, which is just the Christoffel symbol 
standard Einstein field equations in terms of g'^^. The first order Lagrangian reduces to 
the second order Lagrangian by substitution of the solution F^^^ = into the first order 
Lagrangian. 

A canonical analysis of the first order form of the EH action was given for the first 



time by Arnowitt, Deser and Misner (ADM) 2J, |25|. (They also refer to some preliminar 



24 



unpublished steps based on the first order action made by Schwinger.) However, in 
the Dirac procedure was not used and preliminary Lagrangian reduction was performed to 
obtain a reduced Lagrangian with fewer fields than are used in a canonical formulation. To 
do this reduction, the time independent equations of motion are solved to eliminate certain 
fields. 

Straightforward application of the Dirac procedure in the case of the first order for- 
mulation of gauge theories such as Maxwell theory is well-known. In this approach, 
conjugate momenta to all independent variables are introduced and this immediately pro- 
duces an equivalent number of primary constraints as all velocities enter the Lagrangian 



tried to use the Eddington, pure affine, formulation |2l |). so. the symmetry of g"^ and V^^ in a/3 was 
also lifted. (For further developments along this line see |22|.) Einstein considered this formulation as the 
best starting point for possible generalizations of GR. 
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only linearly. From this point, we follow the standard path by considering the conservation 
of constraints in time which produce secondary and higher constraints, until it is possible 
to have all constraints conserved. The AD Maxwell Lagrangian gives 14 constraints, two 
of which are first class and the twelve remaining ones are second class [26]. However, this 
is only a demonstration of the consistency of the Dirac procedure. The next step is the 
elimination of all second class constraints by passing from Poisson brackets (PB) to Dirac 
brackets. The Dirac reduction of Maxwell Lagrangian in first order form is performed in 
Appendix B providing a proof of the equivalence of using the second order and first order 
actions in a canonical analysis. 

A brief discussion of applying the Dirac approach to can be found in j^l where ex- 
pressions for the primary constraints are explicitly given, emphasizing that the first order 
formulation of the EH action in AD results in 50 primary constraints, serving as an illustra- 
tion of the complexity of the Dirac procedure. The number of independent field components 
of g°'^, F^^ in is [d + 1) md dimensions, and introducing conjugate momenta doubles 
the number of phase-space variables. This large number is a way of showing the complexity 
of the EH action, but this is not a real problem, as in the Hamiltonian analysis we separate 
only spatial and temporal indices of fields so that, in this case, we have only nine distinct 
fields for all d: (^^o, (^^fc^ ^fcm^ po^^ ^0^^ ^0^^ T^q, F^, T^- This is not greatly different from 
using four fields in the first order formulation of electrodynamics Aq, A^, Fo^, Fkm- 

In the 2D limit, the first order action is not equivalent to the second order action, which 
is a total divergence . This was analyzed at the level of the Lagrangian in The first 
order Lagrangian in 2D is not a total divergence and its canonical form can be discussed 
just like any other model of 2D gravity j28|]. Moreover, the first order formulation as a 
general field-theoretical construction should be valid in all dimensions, with possibly special 
behaviour in some particular dimensions, but also with some similarities in all dimensions^. 

These considerations have motivated us to perform a canonical analysis of the 2D EH 
action using the Dirac procedure without any a priori assumptions or restrictions. In par- 



^ The possibility of similarities of the 2D limit of the first order form of the EH action with the higher 
dimensional form has to be stronger than is possible in the case of electrodynamics. In the 2D limit of 
the first order form of the EH action we have nine distinct types of fields just as in all higher dimensions 
while in electodynamics there are only three fields in 2D as opposed to four in higher dimensions (as in 
2D, Fkm = 0). 
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ticular, it is important to find the algebra of constraints. In Dirac's analysis of GR |7| the 



PB algebra of constraints 



a 



a 



is non-local with field dependent structure constants 

{Ha (x) , Hb {x')} = Hb (x) S^a {x, x') ~ (ax ^ hx) , (4) 

{Ha{x),H{x')} = H{x)6,a{x,x'), (5) 

{H (x) , H {x')} = h"'' (x) Ha (x) 6,b {x, x') - {x ^ x') . (6) 



This type of algebra (sometimes called a hypersurface deformation algebra [29]) is not 
encounted in 'ordinary' gauge theories. This is not a true Lie algebra, this being the main 
obstacle to canonically quantizing GR. The question that arises is whether this is an intrinsic 
property of GR or the result of assumptions made in the course of analyzing the reduced 
Lagrangian in the approach of Dirac p, 0] and ADM 0, |2^. It turns out that in the 
2D case, the Dirac procedure gives a local algebra of constraints with field independent 
structure constants [30|. In order to preserve 'ordinary' properties beyond locality of the PB 
and to have also off-shell closure of the PB algebra of generators and off-shell invariance of 
;he Lagrangian, we have made a simple linear transformation of the affine connections. In 



30[ , such transformations were expressed in component form but in fact they can be recast 
in the covariant form 

This covariant change of variables is quite different from the usual non-covariant change 
and it provides an alternative covariant formulation of the first order form of the EH action 
which is more suitable for canonical analysis than the form of Q. We have not been able 
to find any particular geometrical significance of the variables but it appears that they 
reflect the dynamical properties of fields of the first order EH action which is richer than the 
geometrodynamics of space-like surfaces. According to Hawking, using a family of space- 
like surfaces is in contradiction to the whole spirit of General Ilelativity and restricts the 
topology of spacetime [s^. (This echoes Dirac's conclusion in ^j, partially cited above.) 
This restriction, imposed by the slicing of spacetime, must be lifted at the quantum level 
I29I 1 ; avoiding it at the outset seems to be the most natural cure of this problem. The idea of 
slicing spacetime originated in the attempt "to recover the old comforts of a Hamiltonian- 
like scheme: a system of hypersurfaces stacked in a well defined way in spacetime, with the 
system of dynamical variables distributed over these hypersurfaces and developing uniquely 
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from one hypersurface to another" 32|- This, although 'reasonable' from the point of view of 



classical Laplacian determinism, is hard to justify from the standpoint of General Relativity 



33 1 . In GR, an entire spatial slice can only be seen by an observer in the infinite future 
I and an observer at any point of a space-like surface cannot have information about the 
rest of a surface. (This actually follows from just the basic principles of relativity, those of 
locality and the finite speed of signals (e.g. see p. 7of ^).) It would be unphysical to build 
any formalism by basing it on the development in time of data that can be available only 
in the infinite future and to try to fit GR into a scheme of classical determinism and non- 
relativistic Quantum Mechanics with its notion of a wave-function defined on a space-like 
slice. This idea also contradicts the canonical treatment of local relativistic field theories 
which do not make any explicit references to the ambient space-time by making use of a 
particular coordinate system or class of coordinate systems 

The change of variables of (|7j) can be used in any dimension and it is quite natural 
to explore this change in higher dimensions with hope that, as in the 2D case, it leads 
to an algebra of constraints that has the form of a Lie algebra or to see how the non- 
locality associated with the "hypersurface deformation algebra" appears in higher dimensions 
without imposing it from outset by choosing a particular slicing of spacetime. 

In the next section we consider the effect of using in place of F^^ in any dimension 
and demonstrate that straightforward application of the Dirac procedure is considerably 
simplified by this choice of variables. After a few simple steps using Dirac reduction to 
eliminate second class constraints, we face sharp discrepancies with previous results [2^ 
obtained by using Lagrangian reduction in which time independent equations of motion are 
used to eliminate some variables. The source of this difference and the conditions under 
which the two approaches are equivalent are analyzed. The next two sections provide the 
full canonical analysis of a simple model both using the Dirac approach (Sec. 3) and using 
the Lagrangian reduction (Sec. 4) in a way similar to in order to illustrate the general 
considerations of Sec. 2. The results are summarized in a conclusion. In Appendix A an 
alternative first order formulation of the EH action in which the variables are used is 
demonstrated to be equivalent to the second order form of the EH action. In Appendix B 



^ The condition that a space-Uke surface remains space-hke obviously imposes restriction on possible coor- 
dinate transformations, thereby destroying four- symmetry. 



8 



we perform the Hamiltonian (Dirac) reduction with the first order formulation of Maxwell 
electrodynamics by eliminating those secondary constraints that are of a special form (this 
is an illustration of what was done in the EH action in Sec. 2 and in a simple model in Sees. 3 
and 4) and prove in this way that, starting from the first order form, one can obtain all of 
the standard results usually derived using the second order form of the action. In Appendix 
C, Lagrangian reduction of the first order form of the EH action in any dimension based on 
the variables is performed in a way consistent with the Dirac analysis. 

II. CANONICAL ANALYSIS OF FIRST ORDER FORM OF THE EH ACTION 
IN ANY DIMENSION 

In this section we discuss the Hamiltonian formulation of the EH action using the gen- 
eralization of the transformation of (|7j) that produces canonical results similar to those of 
ordinary gauge theories in the 2D limit of the first order form of the EH action. The inverse 
transformation of (|7I) in d dimensions is given by 

r^^ = ^ap-^ (^.a + ^^a) (8) 

which upon substitution into Q gives 

L, {g, = {ilp^, - iUh + ^^^A^^.) , (9) 

an alternative first order form of the EH action. This is because the linear transformation 
used for the field F^^ appears in (jS)) at most bilinearly and only linearly in their derivatives. 
It is also possible to prove the equivalence of the first order form © with second order form 
(H)) by solving the equation of motion for and substituting the resulting expression for 
into the equation of motion for g^^^ ^ . As a result, we obtain the Einstein field equations 
without any reference to the affine connection. Actually, solving the equation of motion for 
is simpler than solving that of F^^. (Details are given in Appendix A.) 
However, the main advantage of is that it is extremely well suited for applying the 
canonical procedure. There is now nice separation of components of into those which 
are dynamical and non-dynamical, as the only term with derivatives is of the form 

^ This is similar to what was done to prove the equivalence of l|3Jl to Q in [iof. 
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(Latin indices are spatial and a dot represents a time derivative.) Following Dirac, the first 
step is to introduce momenta conjugate to all fields 

vr.,(^?"^),n^^(e°,),nf (ey . (11) 

Using ©, we immediately obtain the primary constraints 

t:^P ^ 0, nf ^ 0, lif - ^gg^^ ^ (12) 

which equals the number of fields in the Lagrangian. 
The total Hamiltonian is 



al3 
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He = -h-^ - iUh + TT^^^A^^J , (13) 



1 

rf- 1 

where A"^ and A^^ are Lagrange multipliers associated with the primary constraints. 
If the d{d + 1) by d{d + 1) matrix 



Md=({0,0}) (14) 



built from the non-zero PB among the primary constraints (0, G ij^aii, ^^"^ — \f~gg"'"^ 
is invertible, we have a subset of constraints which are second class. Moreover, all these 
constraints are of a special form involving the canonical pair (f?"^, tTq,^) for which n^p ~ 
(see Dirac j^, Appendix B, and for more detailed and general discussion j3|). For such 
constraints, if det Md ^ 0, we can set all momenta TTap to zero and then solve Hq'^ = \/—gg"''^ 
for (7"^ = g""^ ij^") and use this equality to eliminate (7"^ in both the Hamiltonian and the 
remaining constraints. This is Hamiltonian (Dirac) reduction in its simplest form. Dirac 
brackets are equal to PB for all the remaining variables. The use of this reduction is 
shown in Appendix B to lead to equivalence of the first and second order formulations for 
electrodynamics at the level of the Hamiltonian. 

Actually, for Ld it is not even necessary to solve the equations Hq"^ = \/—gg"'" for g'^^ as 
they enter the Hamiltonian in the particular combinations which are present in the second 
class primary constraints and the solution for such combinations are, of course, obvious if 
the condition det Md 7^ is fulfilled. In this case, the canonical analysis of the reduced 
Hamiltonian leads to the form of the gauge transformation of Hg^ and so using the strong 
equality, Hq^ = yf—gg"^, we can immediately find the gauge transformation of g"'^. 
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In 2D (and only in 2D) the matrix ()14|) is singular. The rank of M2 is four and its 
dimension is 6 by 6, so only two pairs of constraints constitute a subset of second class 
constraints which are of a special form meaning they can be eliminated. (For more details 
see Q-) 

Let us denote the number of independent components of a field $ by [$]. At this stage, a 
reduction of the [g""^] = ^d{d+l) fields has been performed by eliminating the canonical pairs 
{g°'^, tTq,^) using the primary second class constraints. The remaining primary constraints 
will produce secondary constraints {xt^) 

nf = {nf.ft} = -^.^ = xf (15) 

with [xf] = ^id~l)d{d+l). 

Explicitly separating time and space indices, we obtain three secondary constraints 

(xr,xr,xr) 

xr = -hT - h'^'Ck - /^""C + (h'-'-CxS]: + h'^^'CxST) , (16) 



^ = xr = -hT - h'Xk - h'Xk + T^/^'^^eiA^r, (17) 

^^ = XT = -/^T - ^h'Vok - 2/^°-e.- (18) 

From the point of view of the Dirac procedure, these three constraints are quite different. 
The distinction that arises between these constraints is not taken into account when they 

agrangian equations of motion in the ADM approach to 



25 



4l|. The matrix of PB of x™" with the corresponding 



are treated as time independent 
the first order Lagrangian of Q 
primary constraints 11™" is non-singular. These constraints form a second class subset and 
this is a subset of the same special form as the part of second class primary constraints 
that have been already considered. Following Dirac reduction, we have 11™" = and C,mn — 
^mn(^a/3) ^Op) ' solutious of the sccoud class constraints xT" = that are now substituted into 
the Hamiltonian and remaining constraints. For the second equation (fT7j) . the matrix of PB 
of xT "with the correspondent primary constraints (H^™) is singular and this subset is not 
purely second class. According to Dirac, we have to find the maximum possible number of 
first class combinations for this subset and only the remaining constraints which are second 
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class can be eliminated. Among the constraints Xfe"^ of ()17p. only one first class combination 
exists and it is xT 

Xf = -h%' - h-'^Ck + h'Too- (19) 

Only the fields are present in (|19|) and they have a vanishing PB with the primary 
constraints H^^. The remaining constraints are again of the special form, so we can further 
reduce our system by eliminating [x^™] — 1 = (rf — 1)^ — 1 fields. Elimination of these fields 
without destroying tensorial notation is performed in Appendix C. 

The last constraint (fTH|) is first class as there are no components of appearing in (fTHjl . 
but only these components give a non-zero PB with the primary constraints H^'^. Dirac 
reduction using pfilll8j) leads to the following number of fields in the reduced Hamiltonian 

[g'''] + la,] - [xD - [xM + 1 = + 2). (20) 

Taking into account this reduction using the second class subset of primary constraints, 
we have only d primary constraints ([n^°] + 1) and d secondary constraints left. 
The secondary constraints are now 

,^00 _ ttOO oTTOOfO oTTOmfO 
Xk — '-'-0,k '^'-h SOA: ^■'■■'■O Smfc5 

xf = -K'k - nreL + KToo (21) 

where we have used the strong equality h"'^ = Uq^ . These are a d dimensional generalization 
of two of the three constraints found in 2D js^, having a simple local PB algebra 

{xT (x) , (y)} = X': (x) S'~' (x-y), {xT, xT} = 0, {xT, = (22) 

and zero PB with primary constraints. At this stage the primary and secondary constraints 
form a first class system. We now continue the Dirac procedure; it leads to the existence of, 
at least, tertiary constraints. 

Already after the first steps of Dirac reduction using ^^/j ^^id g"'^ as independent vari- 
ables, we see that primary and secondary constraints having a local PB algebra with field 
independent structure constants arise and that tertiary constraints must be present in the 
Hamiltonian, which is no longer a linear combination of secondary constraints as in the 
2D case 0J^- This result is quite unlike the previous treatment of the first order EH 
Lagrangian where after Lagrangian reduction the Hamiltonian is a linear combination 
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of secondary constraints with a non-local hypersurface deformation algebra of constraints 
with field dependent structure constants. This has been viewed as an inconsistency in the 



37 



or as an 



constraint algebra and the main obstacle to canonically quantizing GR 
indication of the non-locality of Nature and an inspiration for new ideas, such as promoting 
this algebra to bein g a first principle, more fundamental than the action principle or the 
equations of motion 40|. 

Recently, Kummer and Schiitz have reconsidered the first order formulation of AD GR 
using Cartan variables 0]. Their analysis is based on avoiding the ADM decomposition 
that has been almost exclusively used when discussing tetrad gravity. Their approach also 
leads to tertiary constraints and a local algebra of constraints. 

Before continuing with the Dirac procedure it is necessary to understand why the two 
approaches, Dirac and Lagrangian reduction, that are supposed to be equivalent lead to 
different results. We attempt to answer this question in the rest of this paper. 

First of all, let us note that the presence of tertiary constraints is not in contradiction 
with the number of degrees of freedom. For example, if tertiary constraints are all first 
class and the Dirac procedure is closed at this stage, we have 3d first class constraints. The 
number of fields in the reduced Hamiltonian is d{d + 2) (see (HO))) minus ^d{d + 1) because 
of the first reduction using the second class primary constraints. The result is ^d{d — 3), the 
number of degrees of freedom associated with a symmetric tensor gauge field in d dimensions. 
This expression works only for d > 2; 2D is a special case which cannot be described by 
this relation because there are no second class constraints among the secondary constraints. 
(See for a full discussion of the 2D EH action with g"'^ being treated as an independent 
field.) Of course, having 3d first class constraints is not the only possibility that results 
in the expected number of degrees of freedom, but this demonstrates that the presence of 
tertiary constraints is not inconsistent. Moreover, it is necessary to have tertiary constraints, 
because without solving the first class constraints we have additional independent variables 
and extra constraints are needed to reduce the number of degrees of freedom to the expected 
value. 

Secondly, Lagrangian and not Dirac reduction was used in_|25j. (For a very clear expo- 
sition of this reduction see Appendix A of the review article |4l|.) 

In the Dirac approach, after identifying all second class constraints of the special form 
among the constraints (jKiHlSj) and then eliminating the corresponding variables, we reduce 
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the number of fields (see ()20p) to just d{d + 2) which is 24 if c? = 4. By way of contrast, 



m 



the solution of the 30 Lagrangian equations of motion that do not involve time 
derivatives of any component of the affine connection leads to a reduced Lagrangian with 
only 16 independent variables (see eq.(4.1) of ^25] and also eq.(A.27) of 41]), so that when 
using Lagrangian reduction 34 variables have disappeared after solving only 30 equations 
of motion. This is clear indication that secondary first class constraints have been solved 
in this approach. When using the variables ^2/3? if we were to use the solutions of the first 
class constraints of ()2H) to eliminate fields, we eliminate more variables from the Lagrangian 
than equations of motion that have been solved. This contradicts the Dirac prescription 
for treating constrained dynamical systems and illustrates the importance of classifying 
constraints into first and second class, though the importance of this classification has been 
deemphasized in 42 1. 



The Lagrangian reduction of the Maxwell action written in first order form is, by way of 
contrast, fully justified as all time independent equations of motion in this case correspond 
to second class constraints ^. The equations of motion which do not have time derivatives 
(the Lagrangian constraints) for the first order formulation of electrodynamics are 



km 



§pkm 

giving Fkm immediately in terms of Am- This is in agreement with Dirac reduction, as the 
primary constraints n'^™ ~ (where II'^'" are the momenta conjugate to Fkm) give a subset 
of the second class constraints of the special form which allows for Dirac elimination or 
ensures that the reduced Lagrangian is equivalent to the initial one (see Appendix B). 

Full correspondence between two reduction procedures exists only if we eliminate variables 
using Lagrangian equations similar to (j^Hj) where the field being eliminated is the same as 
the field being varied. (This also is the situation for auxiliary fields in supersymmetric 
models.) Only in this case is the reduced Lagrangian equivalent to the original. We see 
therefore that elimination of variables in the original Lagrangian by merely solving the 
Lagrange constraints may not always be correct. For example, suppose we have an action 
functional S{Q,q) and that the equation 6S/6Q = can be solved for the Q's so that 
Q = Q{q)- This is then substituted back into S and the new action S'{q) = S{q,Q{q)) 



^ The usual references (e.g., |2j|) concerning the similarity of this reduction to reduction of the first order 
form of the EH action are not entirely correct. 
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implies a dynamical equation SS'/Sq = which is the same dynamical equation for q that 
follows from 6S/6q = 0. However, if 6S/6Q = is solved for q instead of Q one does not, 
in general, obtain the same dynamical equations from the action obtained by substitution 
of this solution into the original action. (We illustrate this in Section 4.) 

When using the variables F, it is difficult to compare results of these two approaches 
since the straightforward Dirac procedure is not easy to apply because of the way the La- 
grangian ^ depends on derivatives of F. The subset of second class primary constraints 
used in |^], where the variables F are employed, leads to elimination of some of the vari- 
ables, which affect the rest of the primary constraints and some primary constraints become 
combinations of momenta and not just simply FI^^, as in the case when using the variables ^. 
Consequently not all secondary constraints have the same simple special form structure with 
primary constraints. In Lagrangian reduction, solving those equations of motion without 
time derivatives for auxiliary fields such as is done in ()23p. is not as easy to analyze because 
the last equality in |T5|l is due to the "diagonal" form of terms with derivatives, which is 
not the case when one uses F instead of ^. 

To illustrate this general discussion we do not need to consider the EH action in d dimen- 
sions (using either F or ,^); we just need a simple model in which the first class constraints 
are present and can be algebraically solved in order to compare the results of the two ap- 
proaches. We consider L2 {h,C,) which is a simple model with first class constraints that can 
be algebraically solved ^. It has stronger connection with (g,^) for d > 2 than L2 {g,0- 



3 



with that of L2{g,^) in 36j). 



(To see this, compare the analysis of L2{h,^) in 

In the next two sections we present the complete canonical analysis of this simplest 
Lagrangian using both approaches. We examine the constraint structure and analyze the 
invariance of the action under the gauge transformation that is implied by the full set of 
first class constraints using the approach of Castellani 



^ The non-equivalence of this model to the second order EH action in 2D, the difference between the h and 
g formulations in 2D, etc. are irrelevant since we are only comparing different methods of reduction. 
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III. CANONICAL ANALYSIS OF A SIMPLE MODEL USING DIRAC REDUC- 
TION 



In this section, we give the full Dirac analysis of a slightly modified form of Q which 



has been advocated by Faddeev 



4l|. It allows us to demonstrate the effect of neglecting the 



contributions of surface terms. We consider 

that in the 2D case results in 

Z; = -h^^er, - 2h'ro, - h'Voo - (25) 

where 

+ 2eo\ {hf + h'r^^ - h'Too) + Co^o {h'? + 2/^°^^?! + 2/.°°^) • (26) 

Introducing conjugate momenta to all variables, tTq,/? (^h°'^^ and 11"^ {^Iib) 5 "we immediately 
obtain the primary constraints 

= n^(3 + - 0, Uf ^ (27) 

and the total Hamiltonian 

Ht = H, + A"^<I>„;3 + Al/^n^^. (28) 
Among the primary constraints, we have a subset which is second class as 

{$.,,nr} = A::^. (29) 

We are using the standard fundamental PB for independent fields 

{/."^ TT,,} = A^^, nr} = 5>r„ 

where A"f; = 

I (5° + The constraints ^ are of a special form and, according 

Dirac reduction, such constraints can be eliminated without affecting the PB of all the 
remaining variables. We have now two strong equalities 

no'' = 0, = -^ap (30) 

16 



and as a result, we have the reduced total Hamiltonian 



H^^ = -e^^xV - 2il^xf ~ il^xf (31) 



where using (IHUj) 



Xf = - (/^T - 2/.°Vn - 2/.°°7roi) • (32) 
Conservation of the primary constraints in time leads to the secondary constraints 

n}i = {ni\ h) = Yif = {n?\ h) = , n?" = {n^, h) = xf. (33) 

All secondary constraints have zero PB with the primary constraints and among them- 
selves have the following PB ^ 

[xf, Xf] = xT, Xl'} = -X\\ {xl\ xT} = 2xT- (34) 

The Hamiltonian ()31|) is a linear combination of secondary constraints and, because of 
the PB of the Dirac canonical procedure is completed by the presence of six first class 
constraints for the six remaining canonical pairs (already three pairs Ilg'^ have been 
eliminated) resulting in there being zero degrees of freedom. The secondary constraints 
(j32|) in Lagrangian language correspond to equations of motion obtained by varying such 
non-dynamical variables as those of (fTHj) . These equations cannot be solved for the fields 
being varied and so they are not auxiliary fields. These equations correspond to first class 
constraints in Dirac language and are different in this respect from the algebraic constraints 
arising in Maxwell electrodynamics ()23|). (The effect of performing a reduction by using the 
solution of first class constraints will be considered in the next section.) 

In the Dirac procedure we have to first find all constraints, then eliminate the second 
class constraints. Only after these steps have been performed can we discuss gauge fixing, 
etc. 



® Terms with derivatives must be carefully treated ta king int o account the distributional character of PB 
in the infinite dimensional case, that is, field theory 



SH- 



IT 



To find the full gauge invariance of the action using the Castellani procedure j43|, it is 
important to determine the complete set of first class constraints. If some of the first class 
constraints are solved, they will not be present in the gauge generator and the some of gauge 
symmetries cannot be restored. We can obtain, at most, partial gauge symmetries (e.g., only 
the spatial diffeomorphism) or even possibly the wrong gauge symmetries. 

_43^, is found by first 



The generator G of gauge transformation, following Castellani 
setting G"]^) = Cp for the primary constraints [Cp = ,11^^ and then determining 
^(0) (^) ~ ~ {C'p> ^c} (x) + J dy a'^ (x, y) Cp (y) where the functions a'^ (x, y) are found by 
requiring that |G"q), ifc} = 0. The full generator of gauge transformation is given by 

G (£^ e") = Jdx (e" {x) Gfo) [x) + (x) Gf^ (x)) . 

In our case this leads to the following expression, using the three primary and three secondary 
first class constraints, 

Gis) = l dx [s (-x?^ - ^o'on?" + ^i\n}i) + en?^ 

+^1 {-XV - 24ni^ - 2eooIlf) + e^Ul' + {-xT + '^il'K + '^il^T) + ^^n^] • (35) 

The PB of generators (jH^I) have a closed off-shell algebra similar to that of ordinary gauge 
theories: 

{G{e),G{7])} = G (r^ = C'^^e^rf) ) (36) 

where e'^ = {e^{e),e^{6i),e^{e^)) and the only non-zero structure constants C"^"'' are C^^^ = 
2 = = 1 = _C221^ 6*331 _ ^ _ _^3i3_ ^^ggg ^gflg^^ ^^le structure of the algebra 

of the PB among the first class constraints. More explicitly, these relations are 

r = 2e^r]i — 2£ir]^, ti = er]i — EiT], = e^r] — erf^. 

Now we can find the transformations for all fields appearing in the initial Lagrangian 
that follow from 5 {field) = {field, G}: 

^^u = ^1 - 2£:i^oi + ^^11) 

Seoi = ls-e,eoo + ^'en, (37) 
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Sh^^ = e^h''-e'h''', (38) 
5nu = Ei^i - 2£:i7roi + etch, 

Shqi = l-e^i - ^iTToo + £^11, (39) 



2 

From (jH^ and using the strong equalities (jHUj) . we obtain 



^en = -^i,i-2£ieo°i+^e?: 
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One can easily check the gauge invariance of the Lagrangian of using the trans- 
formations of dSZEHEOl)- It is 

6L', = [h'^e,,, + 1 + ^ - + + ^ (41) 

and so is invariant up to total derivatives. However, a variation of the total derivatives 
appearing in results in a contribution that exactly compensates (jHJ. Keeping the 
initial form L2 of the Lagrangian, we have exact invariance under the transformations of 
()371 1381 HOj) . This illustrates the importance of surface terms in retaining invariance of the 
Lagrangian and shows that the elimination of surface terms can affect its gauge invariance. 
(For a discussion of a similar occurrence in SUSY models, see j^^.) 

The transformations of (p?7| IHH| HlHl can be written in a compact form which is similar to 
one appearing m 

^Cf^ = -^''Up - {ea,Cp. + 4U) , (43) 



^ Due to the transfromations of |4^ we actually were able to recognize the possibility of recasting the 
change of variables (jHJ into covariant form. These were found initially in component form in [30|. A 
similar change of variables is used for the antisymmetric part of the affine connection in generalized 
GR models formulated along the line of [iof. 
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where e"'' is the antisymmetric tensor (e"^ = 1) and is a symmetric tensor with compo- 
nents Coo = £^^5 Cii = ^1) Coi = h^- From using ((71 |H)) we obtain the transformation of 



a/3 



(44) 



— e 



per 



A Ia a ^A ^A 

'^apCpcr — -^^p^avCpa ~ ^oJ^ jipQucr + pyCpa + 2 '^«■'^/3l'Cpcr 



+ (a ^ /5) . 

The Einstein form of the Lagrangian L2 (/?.,r) is invariant under the transfromations of 
and (gH). 

All the usual canonical properties of non-Abelian gauge theories are present in the ap- 
proach outlined here. We have local PB with field independent structure constants, a closed 
off-shell algebra of generators, and exact invariance of the Lagrangian under gauge transfor- 
mations of the original fields. In js^, only the surface term — [h'^^^ap) ^ "^^s added to L2 
and so played role of a generalized coordinate 



IV. LAGRANGIAN REDUCTION BASED ON SOLUTIONS OF FIRST CLASS 
CONSTRAINTS 



Again, starting with the same Lagrangian (j25jl . we find the equations of motion associated 
with the non-dynamical fields to obtain the Lagrangian constraints. These are identical to 
the PB of the primary constraints with the Hamiltonian . We find that 

5L' 



- h]^ - 2h'roi - 2h'roo = 0, (45) 



sen 



§t = + h^r^, - h'Voo = 0, (46) 

si' 



These are the only three equations of motion out of nine in total that have no time 
derivatives. 



This, of course, cannot affect the result as the roles of coordinates and momenta are interchangable in 
the Hamiltonian formulation. The transformations in both cases are the same and only from a purely 
computational point of view might some preference exist. 
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In the Dirac approach, we cannot solve them to ehminate any variables as they are all 
first class constraints. From the Lagrangian point of view, they are not equations for the 
auxiliary fields, since the fields being varied do not appear on right hand side of (j45ll47|) . If 
one were to use them anyway, as is done in ADM approach jj^j, to solve for two variables 
algebraically, such as from and (P7j) 



^00 = ^ {hf + h'r^^) , (48) 

eo, = (^/^T + h^ru) , (49) 

and then substitute these equations back into the original Lagrangian to eliminate ^qq 
and one finally obtains the reduced Lagrangian 

= -h^ru + 2h'''^&, - - (fty - 2hf^ + ft»«^) + s (50) 

where S consists of two total derivatives 

S = In ^ - In ^ . (51) 



(This is similar to eq.(A.26) of 41j.) Note that only two out of the three equations 
can be solved because after substitution of ()48f49|l into pHjl all the variables dissapear, 
converting equation ()45p into a differential constraint that is retained in the reduced La- 
grangian ()5U|). After this reduction, we have five fields out of the original nine left in Zg*^^"* in 
(j3n|) so that two equations of motion have served to eliminate four dynamical fields. This is 
a familiar feature of the Lagrangian reduction in 4D .25]; solving one first class constraint 
leads to the disappearance of two variables in the reduced Lagrangian. In 2^ , after solving 
30 equations, 34 variables disappear from the reduced Lagrangian (as four equations out of 
the 30 that have no time derivatives are first class constraints). 

Moreover, if we perform variation of the reduced Lagrangian ()50|) with respect to ^J'^, we 
obtain 



^--'"'■h-2/.°'^-A"^. (52) 

The variation of the original Lagrangian ()25p with respect to the same variable gives 

^ = -h^' - 2h'%, - (53) 

The use of ()48l49p or any of the other equations of motion cannot account for the difference 
between f|S^ and (|^. Thus, the reduced Lagrangian L2^^ is not equivalent to the original 
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one L'2. If one, despite of this inconsistency, wants to find the Hamiltonian associated with 
this new Lagrangian (j50|l . the Dirac procedure must now be repeated. 
Introducing momenta conjugate to the five fields 

nl^(ey,nSHey,vr„,(/."^) (54) 

leads to five primary constraints 

nP ^ 0, ^ 0, TTn + ill ~ 0, TToi - —ill ~ 0, TToo + -^^u ~ 0. (55) 

Among them we have a pair which are second class of a special form (the second and 
third constraints of ((Hni)^^) that allows us to set IIq^ = and = — tth. Substitution of 
these equalities into the Hamiltonian and the remaining constraints gives the reduced total 
Hamiltonian 

= H^^ + Li + ^^nu) + X'' Lo - ^,^ii) + A}^n}, (56) 



where 



with 



ff^" = -CuxY (57) 



^" = -'^" + 2/,»'^-/,«»^. (58) 

Continuing with the Dirac procedure, we find the secondary constraints. The PB among 
all primary constraints are zero, the only non-obvious one being |7roi + ^tch, ttoo — pn-TTn}. 
The only secondary constraint that arises is 

tll' = {u\\HP} = xl\ (59) 

showing that the Hamiltonian is a constraint. It is straightforward to show that the PB 
among all constraints, both primary and secondary, are zero so that we have four first class 
constraints for four pairs of canonical variables leaving us with no net degrees of freedom. 
At this point everything looks consistent as there are zero degrees of freedom, a local algebra 
of first class constraints and closure of the Dirac procedure. The Dirac constraint formalism 
applied to the reduced Lagrangian gives consistent results even though it is not equivalent 



We can equally take the second and fourth (or alternatively the second and fifth) constraints without 
affecting the final result. 
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to the original theory. As in the previous section, we can find a gauge transformation 
corresponding to this system of constraints. The resuhing gauge generator is much simpler 
than that of It has though the same number of gauge parameters because the reduced 
total Hamiltonian (j56|) also has three primary first class constraints. We find that 



G(e) 



(60) 



The algebra of this generator is closed even off-shell. This generator leads to the gauge 
transformations of the fields 

5ei\ = £i,5np = o, (61) 

^^"=-^5;I^--^;^>^^"=^-^^^"=-' (62) 

'^^ii = ^7;^^ii + ^1.1 - ^1 /;^^T- (63) 

To check the invariance of the reduced Lagrangian (j3n|) . we also need the transformation 
of which can be easily restored by using the strong equality = — tth so that 

< = ^7^^n + ^M-^i^^T- (64) 
Using (I61I6'2I64|) . the variation of (jSHl) is 



=(^7^^n-^i^j^ + ^, (65) 

where S* is a term with total derivatives. Hence the Lagrangian is invariant only on-shell, 
which is a familiar feature of the ADM approach (the gauge generator in AD has a closed 
algebra only on-shell 4^). In the previous section we were able to determine the transfor- 
mations of all fields appearing in the original Lagrangian. It is not possible to do so now; 
by going back we can only restore the transformations of ^qq and .^q^^ using eqs. PHj) and 
but we cannot do this for and Thus using a Lagrangian reduction which is 

based on employing solutions of the first class constraints leads to the gauge transformations 
for only some of the fields in the original Lagrangian and the Lagrangian is invariant only 
on-shell. This is to be compared with having exact invariance of the Lagrangian and gauge 



Of course, after the Lagrangian reduction based on solution of first class constraints we cannot return to 
the variables T because not all the transformations of the fields ^ can be found. It is not difficult to repeat 
such a reduction directly in the Lagrangian when it is written in terms of F and compare with the results 
of the previous section; they will be also different. 
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transformations for all the fields when one uses the Dirac approach to the Lagrangian before 
making a reduction based on all time independent equations of motion. 

Moreover, if we consider the formulation of the 2D EH action using g^^ instead of h"^ 
as independent variables and perform the Lagrangian reduction using (j46p and (j47j) . the 
Lagrangian vanishes identically. To demonstrate this, we present the reduced Lagrangian 
()50|) in the following form 



Lj = -lft"A™-ft»' 



''°'),„]^«n-(*"''"-*°'*°'),];^«n- («) 

If we consider (^"^ to be the independent variable, then h""^ is just a short form for 
^/—gg°'^ and the particular combination that enters ()66|) under derivatives is 

h^V-hV = -l, (67) 

and so, the Lagrangian ()66|) vanishes identically. However, Dirac analysis in this case leads 
to seven independent first class constraints, five of which are primary, and consequently 
there are five parameters characterizing the group of gauge transformations 



V. CONCLUSION 



The second order form of the EH action (Q) in which the metric is the only independent 
dynamical field is invariant under a general coordinate transformations if all terms, including 
the terms with second order derivatives are present. It is possible to apply the standard 
Dirac canonical analysis and to keep simultaneously the effect of all terms when using an 
equivalent first order formulation. The oldest first order formulation which is the closest 
in form to the second order EH action is the affine- metric formulation of Einstein |2ol |. 
This formulation when treated using the standard methods of quantum field theory should 
automatically retain the classical limit. Demonstration that such a limit exists in models 
based on new ideas constitutes a considerable problem in itself ^|. 

The advantage of using this first order Einstein formulation of the action was recognized 
a long time ago by ADM and was used by them as a starting point in their canonical analysis 



of GR [2j, |25[. However, they did not apply the straightforward Dirac analysis and t he) 
performed a preliminary Lagrangian reduction using solutions of first class constraints 



fsee also 



41[). The reduced Lagrangian found in this way is not equivalent to the original 
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Lagrangian and, to quote j47|, "does not represent the full statement of general relativity", 
n the concluding remarks to the last paper of the ADM series (see also remarks in 



48[), the authors suggested that in view of the many ambiguities that could arise in an 



attempt to quantize consistently at reduced level, it would seem more fruitful to return to 
the original Lagrangian (jH)) and try to repeat reduction to the canonical form within the 
framework of quantum theory. To keep the gauge invariance of the original Lagrangian when 
quantizing, all first class constraints must be preserved and second class constraints can only 
be eliminated when they are of a special form, or used to modify the PB by passing to Dirac 
brackets. 

The importance of preserving all first class constraints in the course of the Dirac quan- 
tization was analysed by Ashtekar and Horowitz [49]. They concluded that the Lagrangian 
reduced-space method is likely to yield an incomplete description of quantum gravity. (See 
also the subsequent discussions in |5C|). 

During the two decades between the last paper of the ADM series and the Ashtekar- 



Horowitz analysis, another result of 



47l | got a lot of attention (e.g., see [23] )• This involves 



a geometrical interpretation of what remains after a Lagrangian reduction of field variables. 
The Hamiltonian obtained from the reduced Lagrangian has been emphasized stressing the 
geometrical significance of the reduced set of variables, leading also to a shift back from 
treating four- dimensional Einstein spacetime to just treating space by itself and from con- 
straint dynamics of the full Hamiltonian to geometro dynamics of a reduced Hamiltonian. 
The reduced Lagrangian is invariant only under spatial coordinate transformations and 
the disappearance of some symmetries is a strong indication of the inequivalence of the two 
approaches. The possibility of a 3D geometrical interpretation of the variables appearing 
in the reduced Lagrangian (corresponding to a particular slicing of spacetime) is a demon- 
stration of the inconsistency of Lagrangian reduction because it contradicts the spirit of GR 
and furthermore introduces a restriction on the topology of spacetime This restriction 
originates in solving a part of the first class constraints. This 'freezing' of symmetries is a 
sort of partial gauge fixing used in the ADM treatment [3] right from outset. However, the 
correct procedure is to fix the gauge only after the constraint analysis is performed Q, [s^ . 
This approach also imposes the coordinate conditions that space-like surfaces remain space- 
like as is explicitly pointed out by Dirac . This condition, which restricts the form of the 
general coordinate tranformations, obviously means abandoning four-dimensional spacetime 
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symmetry as is clearly indicated in the conclusion of 

Slicing and the imposition of coordinate conditions also contradict the canonical proce- 
dure, since for relativistic field theories a fundamental tenet of the canonical formulation is 
to not refer to the ambient spacetime 5^ . Any reference to a surface (or a particular subset 
of surfaces) already contradicts this, and it also implies the implicit introduction of extended 
objects into a local theory. It quite likely leads to non-locality, as models explicitly built 
using extended objects (such as string models) are essentially non-local by construction {54 1. 
The common reference to similarities of the hypersurface deformation algebra (^HHI) to string 
models as a sign of consistency of the constraint algebra of the reduced Hamiltonian is really 
a warning sign that GR, which is a local field theory, has somehow been converted into a 
non-local one. 

The Hamiltonian built from the reduced Lagrangian leads to the well-known non-local 
Dirac constraint algebra, which is difficult to quantize as it is not a true Lie algebra. There 
are also problems associated with defining time, finding physical observables, etc.; and as 
a result, numerous attempts to improve this approach by modifying the choice of variables 
01 , reshuffling of the constraints 0, 1^ , etc. have been made. 

However, once reduction has been performed any change of the reduceo? variables or reshuf- 
fling of the reduced constraints cannot cure these inconsistencies, as they are inherent to the 
framework of the reduced Lagrangian. The only possibility of resolving these problems is to 
not abandon the spirit of GR and the standard canonical procedure. The flrst order formu- 
lation of the EH action treated by the standard methods of constraint dynamics preserves 
all symmetries, as it does in ordinary gauge theories, and all results should be reconsidered 
prior to reduction. This was actually suggested by ADM (see Sec. 7-8.1., "Discussion of 
quantization" ) . 

In Isol Isoj l an analysis of the 2D limit of the affine-metric flrst order formulation was 
performed without any a priori assumptions or restrictions such as those used in Lagrangian 
reduction, and the Dirac procedure was applied to see how some properties of GR that 
make it distinct from ordinary gauge theories might appear. However, it turned out that 
all properties of ordinary gauge theories remain manifestly intact if an alternative first 
order formulation is used based on a change of variables involving linear combinations of 
affine connections. This change of variables is easily generalized to any dimension by (jH)) 
and has been employed in this article. The variables C,af3 of (jH)) provide the alternative 
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first order formulation of Q that considerably simplifies straightforward application of the 
Dirac procedure. (For details, see Appendix A.) The first steps of the Dirac algorithm 
(Sec. 2) give results that are different from the ADM analysis which is based on the reduced 
Lagrangian. The origin of such difference lies in using solutions of equations of motion which 
correspond to first class constraints and the reduced Lagrangian obtained in this way is not 
equivalent to the original Lagrangian even at a classical level. (An example of where the 
Lagrangian reduction is equivalent to the Dirac approach is presented in Appendix C.) The 
first class primary constraints obtained in the Dirac analysis ()118|) (Appendix C) cannot 
be eliminated since they constitute a first class subset with secondary constraints (jSI} and 
any further constraint of higher order (e.g., tertiary constraints, either first or second class) 
cannot affect the first class character of this subset of constraints. Moreover, higher order 
constraints do not involve field variables conjugate to variables appearing in the primary 
constraints ()118|) . i.e. the first class nature of primary constraints cannot be changed by 
occurence of higher order constraints. Consequently the Hamiltonian obtained by the Dirac 
approach cannot be reconciled with the Hamiltonian obtained from the reduced Lagrangian. 
The Dirac separation of constraints into first and second class is not merely a technical trick; 
these two classes of constraints are essentially different, as the first class constraints are an 
indication of the presence of gauge invar iance. The knowledge of the gauge degrees of 
freedom is important when quantizing a model and must be kept in the formalism j56| . 

In Sees. 3 and 4, we demonstrated by considering a simple example that Dirac reduc- 
tion and Lagrangian reduction produce different results if Lagrangian constraints are solved 
without appropriate care. At first glance, the simple model which is treated in two different 
ways produces in both approaches the expected canonical results such as a local algebra of 
constraints with field independent structure constants, a closed off-shell algebra of gauge 
generators and the possibility of finding the gauge transformation of all fields. However, 
whereas the Dirac constraint analysis allows us to determine the gauge transformation of 
all fields from the original Lagrangian and to demonstrate exact gauge invariance of the 
Lagrangian, the Lagrangian reduction (based on eliminating variables by use of solutions of 
first class constraints) does not lead to well defined transformations of all the original fields 
and the reduced Lagrangian is invariant only on-shell. 

Moreover, if in more complicated cases the canonical analysis of the reduced Lagrangian 
leads to an algebra of constraints which is not a true Lie algebra, the problem of quantization 
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arises but this may be a problem of the reduced Lagrangian but it is not necessarily a problem 
of the original Lagrangian. 

In our simple 2D example we obtained quite different results using the two approaches 
showing that in the general case the canonical analysis of what is obtained after Lagrangian 
reduction can lead not only to gauge invariance on-shell but also, for example, to non-locality 
of PB, a closed algebra of generators only on-shell [43 1, and consequently to a wrong or, at 
most, only partially right description of the initial Lagrangian. We thus feel that for the 
EH action in higher dimensions it is natural to expect there to be an even more drastic 
deviations between the two approaches. What are they? We are not going to speculate 
here about all possibilities, but we hope that we have been able to convince the reader 
that the existing canonical formulation of the first order EH action has been obtained in 
a non-canonical way and its reduced, geometrodynamical, formulation is not equivalent to 
the original EH Lagrangian. This gives rise to the very important question of what we are 
trying to quantize in canonical qravity. Is it the full Einstein GR theory or only the spatial, 
geometrodynamical, part of it? 

The solution of the first class constraints in the first order formulation is somehow 
related to a partial neglect of surface terms and imposing coordinate conditions in the second 
order analysis jr^. It is natural to ask about the connection between solving first class 
constraints in first order formalism and breaking relativistic invariance by integrating out 
second order derivatives in the second order formulation. The full answer to this question 
can be given only if a generalization of the Dirac procedure is possible that allows us to 
deal directly with accelerations present in the second order EH action. Qualitatively, we 
exp ect that the term linear in acceleration corresponds to a primary constraint (see 

which initiates a chain of higher order constraints. The elimination of such a term 
corresponds to cutting off the first term of a chain of constraints and has the same effect as 
solving a first class constraint in the first order formulation of the action. 

The numerous problems associated with canonical geometrodynamics are quite likely just 
problems of using the reduced Lagrangian (with its reduced symmetries), not an intrinsic 
characteristic of GR, and can actually be considered as an illustration of the fact that having 
only spatial symmetry is not enough for a consistent formulation of GR. The simplest and 
most natural possibility for resolving these problems has not been fully explored: instead of 
trying to improve the reduced formulation or attempting to find some new physics in the in- 
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consistencies of canonical geometrodynamics, one should try to find a canonical formulation 
of GR by applying the Dirac procedure to its first order formulation without any a priori 
assumptions or restrictions. The use of an alternative first order formulation given in Q 
that is based on a generalization of the transformation found in the 2D limit of the action 
provides an example of how Dirac or Lagrangian reduction can be performed consistently. 
Possibly the use of these new variables is not sufficient to ensure a canonical form of GR that 
allows for quantization, and further modifications are needed in order to find a first order 
formulation that preserves all the properties of ordinary gauge theories in higher dimensions. 
In particular, we want to find a formulation that leads to a local algebra of constraints with 
field independent structure constants (as in 
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36| and eqs.(j22ll21I))- The existence of such 



an algebra is needed to pass a crucial consistency test (for Hamiltonians) recently empha- 
sized in Q]. We believe that all possibilities have to be explored in this direction before any 
new physical hypothesis is introduced and before the question: "spacetime or space?" (once 
answered by Einstein) can be posed again 

We would like to conclude our discussion by the epigraph to the second lecture, "Ge- 
ometrodynamics" , in the course on Canonical Quantization of Gravity at Banff Summer 
School [33]: 

There is only the fight to recover 

what has been lost 
And found and lost again and again 

T.S.Eliot: Four Quartets 
East Coker, 186-7. 
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VII. APPENDIX A 

Here the proof of the equivalence of the first order formulation L(i{g,C,) defined in (jH)) 
with Ld (g) of ((H) is presented for the case d ^ 2. 

The variation of Ld {g,C,) with respect to g^^^ gives the standard result 

SLd ( . — . 1 



bg^"" 

where 



g^f^--^gg'-Pg^^R^^{^) (68) 



From ()68p it immediately follows that 



R^v - \g^.uR = (70) 

or alternatively, by using the inverse of the expression multiplying R^is in 

"^^Z - Iv^g'^'g.u) ' = (a^^: - -^g'^g,. ] (7i) 



(not defined in 2D), we obtain 

R,. (0 = 0. (72) 

R^fj (^), as in the case when considering Ld{g,T), now has to be expressed in terms of 
g"^. Varying Ld with respect to we have 

^ = -Z^;;^ - h'^^e,. - h^%u - ^ {h^^e.xK + ^"''^A^.^) ■ (73) 

This equation is easier to solve than the analogous equation for T . First, we can obtain 
the trace of Multiplying (|73|l by h^p (where haph^'^ = 6X) we obtain 

^A^ = -l^Krh^J- (74) 
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We see that here, as when using the variables T, 2D case is special. Substitution of (f7i)) 
into (fTSj) gives 

h'"'e,u + h^'ilu = Dl" (75) 

where 

D7 = -h"; - ^^^^ {h^^sp + h^ps:) h^^h^. (76) 

Multiplying (f75j) by h^phptj and performing a permutation of the indices a, /3, v we obtain 
three equations; we add two (with the permutations (/5, a, z/) and (z/, /3, a)) and subtract the 
third (a, i^, /5), so after multiplication by we obtain the solution for ^"^^ 

C = \ K,D^^ + /^.M^r - . (77) 

Substitution of (f77|) into (fTOj) or (f?^ gives the Einstein equations for free space. (No 
reference to F has been made.) Similarly, if we substitute this solution into the Lagrangian 
Ld{g,C,) of (jni), we obtain the reduced Lagrangian which is equivalent to the second order 
form of EH action L^i (g) (including terms with second order derivatives). 

The appearance of explicit dimensional dependence in ()74l76p seems to be inconsistent 
with using the Christoffel symbol^^. To resolve this, let us consider the trace of F expressed 
in terms of C,- Using (jH)) we find 

= -^C'a' (78) 
so that upon substitution of into (|7H|l and remembering that is only short for 
\/—gg'^^, we obtain 

= ^h^phf = ^-j=9c., {^99'-')^^ = -\9.,9f (79) 



-9 

which is a well-known expression. Similarly, the general case for arbitrary F^^ can be 
demonstrated using (jH} and (fTTj) . 

VIII. APPENDIX B 

As an illustration of Dirac reduction (used in Sections 2-4) that employs elimination of 
only second class constraints that are of a special form, we prove the equivalence of the first 



The expression for it has the coefficient i in any dimension. 
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and second order formulation of Maxwell electrodynamics at the level of the Hamiltonian. 
The first order form of the Maxwell Lagrangian is 

Lm = {d^A, - d,A,) F'^^ + (80) 

where and F^,^ = —Fy^ are treated as independent fields. This formulation is equivalent 
to the standard second order form. This is obvious at the Lagrangian level, as the auxiliary 
field F^jj can be easily eliminated. It is also not difficult to prove this at the Hamiltonian 
level by reducing the Hamiltonian that corresponds to (|8(jp to the standard one by using the 
Dirac procedure. 

Introducing momenta conjugate to all fields tt'^ (^/x), H^,^ {F^^) we obtain 

6 (doA^) 

and 



In AD, equations (j811 182j) give ten primary constraints 

r = T^^- ~ 0, = H^. ^ (83) 
and the total Hamiltonian is 

Hp = H, + X^r + ^^"^^.u (84) 
where A^, A^'^ are Lagrange multipliers and 

= dkAoF'^ + dkArnF''^ - ^F.oF'^o - (85) 

The proof that the Dirac procedure closes can be found in 2^. The first order formu- 
lation produces 14 constraints (2 first class and 12 second class). However, we will proceed 
differently by eliminating step by step the second class constraints that are of a special form. 

The non-zero fundamental PB are 

{A„ n^} = S;, H,.} = i {6^^S: - S;st) . (86) 
It is not difficult to calculate the PB among primary constraints. The only non-zero 



brackets are 
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so that there is a second class subset of primary constraints because 

{vr'=-F^°,no„} = ^5t (88) 

This subset is of a special form and allows one to solve these constraints leading to the 
reduced Hamiltonian (with a reduced number of variables) after substitution of 

Hofc = 0, F'^o = TT^ (89) 

into the original Hamiltonian as well as into the remaining constraints. After the first stage 
of reduction, we have 

H^;^ = Hi'^ + Ao0° + A'^'^km (90) 

where 

Hi'^ = d.Aon' + duA^F^^ - - \f^^F^^ . (91) 

We have now only seven independent fields and four primary constraints 

00 = vr° ^ 0, <l>fc^ = Hfe^ ^ 0. (92) 

They have zero PB among themselves and conservation of these constraints now has 
to be considered. The conservation of the primary constraints gives rise to the secondary 
constraints x° and Xfcm 

= {7r°,i/«} = {.\ -Ao9.7r^} = 9,vr^ = / (93) 

and 

^km = {^km, H^^^} = ~ (dkA^ - d^Ak) + ^Fkm = Xkm- (94) 

It is obvious that the constraints Xkm constitute a second class subset of a special form 
with and, as was done at the previous stage, can be eliminated by solving them for 11^^ 
and 

Ilfcm = 0, Fkm = {dkAm - dmAk) . (95) 
Upon substitution of (|95|l into H^^\ we obtain 

= HP + Ao0° (96) 

with 

hP = dkAon'' - ^TTfevr^ + ^ (dkAm - d^Ak) {d'^A^ - d^'A'') (97) 
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which is exactly the standard Hamiltonian in the second order formulation with only four 
fields and two first class constraints. This completes the proof of the equivalence between 
the two types of reduction at the pure Hamiltonian level. It is important to note that here 
we never 'solve' first class constraints, unlike what occurs in the EH action when treated 
using the ADM formalism and so the problem associated with ADM reduction do not arise. 

IX. APPENDIX C 

We perform Lagrangian reduction of the action of Q in a way that is consistent with the 
Dirac procedure by eliminating only non-dynamical fields by solving equations of motion 
with respect to fields used in the variation that leads to these equations. The Lagrangian 
after a complete separation of spatial and temporal components, can be presented as a 
sum of terms 

where the purely dynamical part is 

= {h'-eo,eom + 2h''eo,eoo + /^°°eo°oC) ■ (99) 

The term with the non-dynamical field ^qo is 

L2 = h'X,^, - 2^00 {h'Xk + , (100) 

while the term with the non-dynamical field is 

1 2 

and finally the term that is at least linear in non-dynamical field is 

OA; 00 fc 1 00 OA: Ok 

Li = 2h "^^Orn,k ~ h iornCok + ^ _ ^Op^Oq "~ 2/l C,OqC,pk + ^ _ ^Op^kq (102) 

9 9 
Oh^kcP cO I ^ uOkcP fO cy^kmcP cO I ^ hOOcP dO 

SOA;SOp "T d _ I 'SOpSOfc SOmSpfc "r ?0p?00- 

The 2D limit of (jHEI) is obtained by setting d = 2 and putting all spatial indices equal to 
one, giving (j2Sl)- 

We have three "non-dynamical" fields among the fields ^ (i.e., fields that enter the La- 
grangian without any time derivatives): .^qo' '^om ^mn- The first field enters only linearly 
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and cannot be eliminated. (This term ()1U0|) corresponds to a first class constraint in the 
Dirac approach.) Variation of ()102p with respect to gives 

- '^h'Xa + -^h^'CoA = DT {h-^ eL) • (103) 

The left side of this equation is not invertible and not all components can be eliminated 
because the trace of the left side is zero. To preserve the tensorial character of variables 
in the reduced Lagrangian it is better to introduce an extra (pure auxiliary) field 6 by 
performing a change of variables in the following term 

I^"°X?S,-^'>°X''-i^«''. (104) 

Introducing of this extra field 6 allows us to solve ()103p for all components of ^g^. Vari- 
ation of ()102|) . taking into account ()104|) . gives 
1 



b _ 



(105) 

Substitution of this leads to the reduced Lagrangian (with the (rf — 1)^ components of 
,^om completely ehminated) 

= h-^i% + (C) + L« (4; + Li'^ {0; C) + l':^ (eL; C,) (loe) 

where now 



k d-l 

4'^ = L2; (108) 

4'^ = ^ (/^"'"^.rn - ^/^'^'"eL + oh'Xo) ; (109) 

7-(l) ^ LmnAfc _ km(:P (:q , 1 fcmfP eg _ Oh'^^ \ ^^f^^ \ ^ i^Om [ ^ ^ c9 ) 

V / ,fc V / ,m 



^00 SOpSmA; ^ ^OO ^OfeSimij ^ ^ ^qO >>kp^nm \-^-^^J 

9 !, fem L On 9 

^ cO I /^mfctO I LOmtO \ 
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Note, that the terms quadratic in 9 have cancelled out and that the terms linear in 9 lead 
to the additional first class constraint ()109p . The terms quadratic in are multiplied by 

!,Ofc !,0m, 

where e'^'" has the property: e^'^hmn = ^n- can now reduce the Lagrangian further by 
ehminating The variation of (jllUp with respect to gives 

e^'ila + e'%la - -^a, {e'X + e^^5^) = (1 12) 



where 



1 uOb 1,0b 1,0c l,kbl,Oc 

+-^^''eok^i - T^^S f/^r + h'-ek^ - + - c) . 



d-l « - 1 /lOO 

The solution of ()112p is similar to solution of (fTSj) appearing in Appendix A. Multiplying 
fjll2j) by /ifcc, we obtain the trace 

= (114) 
and substituting of ()114|) into ()112|) leads to 

e'\L + e'^^eL = + ^^/^P,^r (e'^'^: + e^^^^ = ^f- (US) 

Multiplying pi5j) by hrbhgc and performing a permutation of the indices r, s, a (as in 
(fTSj) ) and then multiplication by Ae"'^, we obtain the solution for 



2 

or in terms of (|113|) 



C = \ IkA" + h^D^ - e"'h Ji,.M (116) 



(117) 



Substitution of ()117p back into the Lagrangian L^]^^ produces the reduced Lagrangian 
with the non-dynamical fields all absent. We can use this reduced Lagrangian to pass 
to a Hamiltonian formulation which is different from the ADM-reduced formulation, as all 
components of are still present (and not only its spatial components). This should lead, 
in principle, to a restoration of full gauge invariance that will involve all components of g""^ . 
The solution of the equations of motion corresponding to the first class constraints in the 
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Dirac approach leads to a non-equivalent reduced Lagrangian with a loss of the possibility 
of restoring full gauge invariance of the initial Lagrangian. This is illustrated for d = 2 in 
Sec. 4. We note that our elimination of the non-dynamical variables produces an alternative 
formulation of the Einstein-Hilbert action that is linear in time derivatives of the dynamical 
fields and well suited for application of the standard Dirac procedure. 

The Dirac analysis applied to the reduced Lagrangian L^p gives (as in Sec. 2) the second 
class primary constraints 

as well as two first class primary constraints 

n^o^O,7rf«0, (118) 

where n°°, vr are momenta conjugate to fields ^qq, 9. Conservation of the constraints of ()118|) 
in time (using ^\()()\ ll()f?|) ) leads to secondary constraints which are equivalent to (^1]) . 
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